Dynamo action at low magnetic Prandtl numbers: 
mean flow vs. fully turbulent motion 
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We compute numerically the threshold for dynamo action in Taylor-Green swirling flows. Kine- 
matic calculations, for which the flow field is fixed to its time averaged profile, are compared to 
dynamical runs for which both the Navier-Stokes and the induction equations are jointly solved. 
The kinematic instability is found to have two branches, for all explored Reynolds numbers. The 
dynamical dynamo threshold follows these branches: at low Reynolds number it lies within the low 
branch while at high kinetic Reynolds number it is close to the high branch. 

PACS numbers: 47.27.eq,47.65.+a91.25w 



The magnetic field of planets and stars is believed 
to be the result of a dynamo instability originating in 
the flow motions inside their electrically conducting fluid 
core. Dynamo occurs when induction due to motion over- 
comes diffusion corresponding to a threshold in the 
magnetic Reynolds number of the flow (Rm = UL/ij, 
with U and L characteristic velocity and length scales 
of the flow, and 77 the fluid's magnetic diffusivity) . For 
liquid metals (as molten iron in the Earth core, or liq- 
uid sodium in laboratory experiments Q), the kinematic 
viscosity v is several orders of magnitude lower that the 
magnetic diffusivity r\ - the magnetic Prandtl number 
Pm = v/rj \& often of the order of 10 -5 or lower. As a 
result, the kinetic Reynolds number of dynamo generat- 
ing flows tends to be very high. Indeed, the relationship 
Rv = ULjv = Rm/Pm implies that critical values R C M 
of the order of few tens are associated with Reynolds 
numbers in excess of one million. This remark prompts 
two initial questions: (i) can a dynamo instability de- 
velop from a fully turbulent flow?, and (ii) what is the 
evolution of the instability threshold as Ry grows? We 
partially addressed them in a previous numerical work Q 
for a flow generated by a deterministic Taylor-Green forc- 
ing at large scales. Lowering the magnetic Prandtl num- 
ber Pm from 1 to 10~ 2 , we established that the laminar 
R C M value undergoes an eightfold increase as unsteadiness 
and small scale motion develop, and that once turbulence 
is fully established the threshold R C M tends to saturate 
to a constant value. Dynamo action is thus preserved 
for flows having a well defined large scale geometry, a 
question still open for turbulent motion without a mean 
flow Q . Numerical predictions of a dynamo threshold in 
realistic natural or experimental conditions are still out 
of reach; numerical simulations cannot get close to the 
parameter range of planetary bodies or laboratory ex- 
periments (even if one ignores the question of boundary 
conditions, a DNS study at Pm ~ 10~ 5 would require a 
resolution in excess of (10 5 ) 3 grid points). Nonetheless, 
the fluid dynamo experiments in Riga and Karlsruhe 
found the onset to be remarkably close to the values pre- 



dicted from numerical simulations based on the mean 
flow structure and this despite the fact that the cor- 
responding flows are quite turbulent. It has led several 
experimental groups seeking dynamo action in less con- 
straint geometry, eventually leading to richer dynamical 
regimes, to optimize the flow forcing using kinematic sim- 
ulations based on mean flow measurements 0, - with 
the advantage that mean flow profiles can be measured 
in the laboratory. It is thus of interest to test the validity 
of this procedure in numerical experiments. 

In this Letter, we compare numerically the dynamo 
behavior as simulated from the dynamical magnetohy- 
drodynamics (MHD) equations to the result of kinematic 
calculations in which the velocity is fixed to its time aver- 
aged profile. In the fully dynamical problem, we integrate 
pseudospectrally the MHD equations : 

^+vVv = -VP + jxB + vV 2 v + F (1) 
at 
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together with V • v = 0, V • B = and a constant mass 
density is assumed. Here, v stands for the velocity field, 
B the magnetic field, j = (V x B)//z the current density 
and V the pressure. In our case, the forcing term, F, 
responsible for the generation of the flow is chosen to be 
the Taylor-Green vortex (TG) Q: 
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with fco — 1- For a given fluid viscosity and forcing am- 
plitude, we first allow the flow to settle into a statisti- 
cally steady state, in a nonmagnetic phase (Eq. ^ with 
B = 0). We either use Direct Numerical Simulations 
(DNS) for flows at Pm of order one, or Large Eddy Sim- 
ulation (LES) schemes with an effective viscosity v e s [Hj 
for Pm lower than about 0.1. In the results reported here 
(cf. Table I), the forcing flow amplitude is kept equal to 
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1.5, and the Reynolds number is increased by lowering 
the fluid's viscosity. Once the steady state is reached, a 
seed magnetic field with energy 10~ 20 , evenly distributed 
among Fourier modes, is introduced and the MHD equa- 
tions are integrated for several choices of the fluid's mag- 
netic diffusivity. For each, we compute the exponential 
growth rate ob = d(lnEB)/dt. The critical R C M dyn , for 
which as changes its sign, is thus obtained at constant 
Rv in a dynamical process for which the velocity field is 
a true Navier-Stokes flow. 

In the kinematic study, we first compute the mean flow 
as an average in time of the velocity fields generated in 
the dynamical runs: 

U(r) = i|dtv(r,t), (4) 

where T is a time scale chosen much larger than a typical 
eddy turnover time at the flow integral scale Ld yn ' ?nl = 
Udyn/Ld yn (Udyn- characteristic velocity). Note that T 
should exceed the magnetic diffusion time T^j = -RajTnl- 
In practice, an instantaneous velocity field is extracted 
every time interval AT to increment the running average 
of U(r). We have used AT ~ T NL /100, and T > 200T NL . 
To save computer time, the averaging is done during the 
linear growth (or decay) phase - hence in the absence of 
a Lorentz feed-back in the Navier-Stokes equation. The 
induction equation (2) is then solved with v(r, t) kept 
equal to U(r), in search of growing solutions B(r, t) = 
cxp (cr B m t) B(r). We note that the mean flow defined in 
Eq. 01 is not a real flow in the sense that it is no longer 
a solution of the Navier-Stokes equations; for instance, 
it does not have a viscosity. We chose to attribute it 
the viscosity of the generating dynamical run. Then one 
could compute an associated kinematic Reynolds number 
Rv,kin = fkinTkin/^, but we have chosen to represent all 
Reynolds number variation as a function of Rv = Rv,dyv 
(Figs. 1 and 3). Critical magnetic Reynolds numbers 
R C M dyn and R C M kin which are computed from each field's 
characteristic lengths (Fig. 2). 

Before analyzing the dynamo behavior, we first com- 
pare characteristics of the dynamic and time-averaged 
velocity fields. Their spectra are shown in Fig.^a), for 
fields originating in the DNS calculation at Rv,dyn = 670. 
While the dynamical flow has a typical turbulence spec- 
trum, the time-averaged field is sharply peaked at the 
size of the Taylor-Green cell. As the Reynolds number 
varies, the quantities for the average flow remain con- 
stant, while they do vary for the dynamical field. For 
instance, the flow integral scale is shown in Fig.^b). It 
is computed from the kinetic energy spectra 

L _/ JdkEy(k,t)/k \ 

2tt \ JdkE v (k,t) J T 1 ' 

where (-)t stands for averaging in time. For Reynolds 
numbers less than about 100, Ldyn tends to be larger than 
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TABLE I: Parameters of the runs: grid resolution, fluid vis- 
cosity (effective value for the last run which use LES), in- 
tegral length and velocity scales, kinetic Reynolds number 
Rv,dyn = Ldy-nUdyn/v, critical magnetic Reynolds number 
^M,dyn f° r the dynamical problem and for the kinematic one: 
R-M,kin f° r the onset of the the first dynamo mode, Am kin the 
value at which it is no longer a growing solution and -Rju.kim 
the onset value for the second dynamo mode. The kinematic 
magnetic Reynolds numbers are computed using the kine- 
matic integral velocity and length scales, and the magnetic 
diffusivity: -R A /, kin = L^JU^/rf ■ For all runs, we have 
measured Lkin ~ 7r and [/kin ~ 3.0. 




FIG. 1: (a) kinetic energy spectra for TGI (y = 0.01); 
Ev,dyn(k,t — T) (solid line), -Ev,km(fc) (dotted line), (b): inte- 
gral length scales Ldyn and Lkin, normalized by the size of the 
unit TG cell, versus the flow Reynolds numbers Rv = Rv,dyn- 

the size (n) of one Taylor-Green vortex. At higher Ry's, 
the turbulent flow has an integral length scale clearly 
confined within the Taylor-Green cell. On the opposite, 
the mean flow has Lkin ~ t and [/kin ~ 3.0 at all i?y's. 

We now turn to the dynamo generation. In the fully 
dynamical problem, the dyn (Rv) curve - Fig. 2(b) - 
displays an initial increase, corresponding to the devel- 
opment of turbulence, followed by a plateau Q. The 
dynamo threshold is then independent of the fluid's vis- 
cosity. This is understood if we recall that the dynamo 
is governed by lar ge s cale motions : indeed, in TG flows, 
we have reported [llj that in the early stage all Fourier 
modes grow at the same rate, but after a few hundred 
large scale turnover times only the modes near the TG 
forcing keep growing while the small scales arc quenched. 
Besides, it is a common observation in fluid mechan- 
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FIG. 2: (a) Growth rates for the kinematic dynamo generated 
by mean flows computed for decreasing values of the fluid's 
viscosity. The intersections with the <r_g = axis define the 
kin vaules reported in Table I. (b) Evolution of the critical 
magnetic Reynolds numbers R C M with Rv- 



ics |12( that the geometry of a shear flow becomes in- 
dependent of the Reynolds number at high Ry, except 
for viscous sub-layers near boundaries. One thus expects 
that, for flows having a well defined geometry, the dy- 
namo persists at all (small) magnetic Prandtl numbers, 
with a finite value of R C M as Rv — > oo. See also ^| and 
references therein for a recent discussion. 



For the time-averaged How, we have observed the ex- 
istence of two dynamo branches in the kinematic sim- 
ulation of the induction equation - a behavior already 
noted for the ABC flow 0. As shown in Fig. 2(a), 
the kinematic growth rate is positive in the interval 

[^M.kin'^M.kirJ ~ [ 22 > 51 L and then a g ain for #M,kin > 

kin- Beyond kin , the growth rate seems to be a 
monotonously increasing function of Rm- The first dy- 
namo window is found to be essentially independent of 
the kinetic Reynolds number of the dynamic flow that has 
generated the time-averaged field. For the higher branch 
it varies within 15% of a mean value R M km ~ 158 (with 
no systematic trend in the range explored in our study). 
We compare in Fig. 2(b) the evolution of the critical mag- 
netic Reynolds numbers. At low Rv, we observe that the 
dynamo threshold for the dynamical problem lies within 
the low Rm dynamo window for the time-averaged flow. 



For Ry larger than 200 the dynamical dynamo thresh- 
old lies in the immediate vicinity of the upper dynamo 
branch (high R C M mode of the time- averaged flow). 

Let us now compare the structure of the resulting dy- 
namo fields. The magnetic energies and corresponding 
integral length scales are shown in Fig. 3. For the dy- 
namical dynamo, energy is distributed in a broad range 
of scales. For the kinematic dynamos evolved from the 
time averaged flow, we observe that for the low Rm mode 
the energy is strongly peaked at large scales, while it is 
more evenly distributed in the case of the high Rm mode. 
This behavior is also reflected in the evolution of the inte- 
gral magnetic scales, shown in Fig. 3(b). The kinematic 
low mode grows a dynamo essentially at scales larger than 
the Taylor-Green cell {Lb ~ 1.6n) at all Ry's. The high 
mode on the other hand grows within a TG vortex. The 
dynamo selected by the dynamical How switches between 
these two behaviors. At low Ry it grows with an inte- 
gral scale larger than the TG cell and we observe that 
Lb, dyn ~ -^s,kin- At high Ry the magnetic integral scale 
is about half the size of the TG cell. The peak in the 
magnetic energy spectrum at smaller scales in the dy- 
namic runs, as well as the smaller integral scale, suggest 
that turbulent fluctuations can play a role in the dynamo. 
This is in agreement with |lll| , where both large and small 
scales were observed to cooperate in the dynamo process: 
the growth of small scales allows the dynamo to quench 
velocity fluctuations and the magnetic field to grow faster 
at large scales. 




FIG. 3: (a): Magnetic energy spectra for v = 0.01 (DNS 
at Rv~G70); (solid line) dynamical run, (short and long 
dashed lines) low / high Rm kinematic dynamo modes (in 
this regime the spectra are computed as averages during the 
growth phase, normalized by the mean energy), (b) Evolution 
of the integral scales for the magnetic field, computed from 
the spectra. 

The structure of the dynamos can be explored further 
with visualizations of the isosurfaces of the magnetic en- 
ergy. Fig. 4(a) and (c) correspond respectively to the 
low and high kinematic cigenfunctions, while Fig. 4(b) 
and (d) show the dynamical fields, with the magnetic 
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energy rescaled and averaged in time during the linear 
growth phase (i.e. (Em(t, t)/£^M(i))r)- Both a low 
(R v = 76.74) and a high (R v = 670) Reynolds num- 
ber are shown. One observe a very good correspondence 
between the low Ry dynamical mode and the kinematic 
low eigenfunction; indeed, at low Reynolds number the 
flow is laminar with small fluctuations about its mean. 
In this regime the dynamo is mainly generated in the 
shearing regions in between the TG cells 0, ^J. At 
high Ry , we compare the dynamical growing dynamo to 
the structure of the high Rm eigenfunction. In Fig. 4(c) 
we recognize the 'twisted banana shape' structure of the 
neutral mode that underlies the a — Q dynamo of von 
Karman flows (recall that at high Ry the TG flow in 
each cell is very similar to the VK swirling flow) 0, . 
For the dynamical flow at high Ry the TG cells are no 
longer as coupled as they are at low Ry, and in Fig. 4(d) 
one does not observe the clear pattern of the kinematic 
cigenmodc. This can be due to the fact that the mag- 
netic energy is displayed here during its initial growth. 
As seen in the spectrum- Fig. 3(a) - the magnetic energy 
grows at all scales; it is only in the non-linear phase that 
the magnetic energy is eventually dominated by the large 
scales. However, at all -fiV's the dynamo grows predomi- 
nantly in the planes 7r/2 and 37r/2 which cut through the 
center of the Taylor-Green cells. 




FIG. 4: Isosurfaces of the magnetic energy, (a): TGI kine- 
matic low eigenmode; (b): corresponding dynamical run at 
Rv = 76.74; (c): high kinematic eigenmode; (d): correspond- 
ing dynamical run at Ry = 670. In each plot two isovalucs 
have been plotted at 50% and 75% of the maximum magnetic 
energy - except for the high kinematic modes in which the 
levels are 25% and 50%. 

We have observed a similar behavior for flows forced 
with TG at fco = 2 (existence of at least two branches 
for the kinematic problem, with a transition of the dy- 
namical dynamo threshold from the lower to the upper 



branch as turbulence becomes fully developed), but re- 
sults may differ for other types of geometries or flow forc- 
ing . However we propose that the observations made 
here, particularly in the large Ry limit, may have some 
relevance for laboratory experiments. For instance, we 
find that once turbulence is fully developed, the dynamo 
threshold is well approximated by the (high branch) kine- 
matic value computed using the time averaged flow. This 
finding is in agreement with the observations in the Riga 
and Karlsruhe experiments, where the mean flow struc- 
ture was optimized to favor dynamo action at low R\[. 
In addition the threshold reaches a limit value when Ry 
increases, in agreement with kinematic simulations using 
von Karman mean flows measured in the laboratory ||: 
the threshold was observed to be Ry independent for 
Ry exceeding about 10 5 . Hence one may expect that 
kinematic predictions based on hydrodynamic measure- 
ments in laboratory prototypes can be relevant for exper- 
iments. This is also of interest for the numerical study 
of natural dynamos where a fully resolved description 
of the fluid's motion in the correct range of parameters 
is out of reach 01 ■ There are however reasons to be 
cautious. The main concern lies in the observation that 
fully turbulent flows in confined volumes are not station- 
ary. Long time dependence in the large scales velocity 
fluctuations have been observed fl8t Also, the ef- 

fect of turbulent fluctuations on the threshold need fur- 
ther study. While recent works have show that a large 
scale incoherent noise may increase significantly the dy- 
namo threshold [2jj , small scale fluctuations can also be 
a source of dynamo action P, [2l| . 
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